Structural and Electronic Properties of Oxidized Graphene 
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We have systematically investigated the effect of oxidation on the structural and electronic proper- 
ties of graphene based on first-principles calculations. Energetically favorable atomic configurations 
and building blocks are identified, which contain epoxide and hydroxyl groups in close proximity 
with each other. Different arrangements of these units yield an LDA band gap over a range of a 
few eV. These results suggest the possibility of creating and tuning the band gap in graphene by 
varying the oxidation level and the relative amount of epoxide and hydroxyl functional groups on 
the surface. 



The prospects of graphene-based nanoelectronics [J, Q 
have stimulated extensive research activities in recent 
years. Pristine graphene, with two linear bands cross- 
ing at the Dirac point, is a zero-gap material. There- 
fore, much effort has been devoted to creating an energy 
gap in CTaphene-based materials for device applications 
In particular, an energy gap can be 
achieved through either nanopatterning or chemical 
functionalization 0, H, ■ The latter has a greater ad- 
vantage because of the ease to scale up in the production. 
A recent experiment has successfully generated a metal- 
semiconductor-metal junction using epitaxial graphene 
and a single functionalized graphene sheet (FGS) 
Analogous to the combination of Si and Si02 in the cur- 
rent generation of microelectronics, the FGS has the po- 
tential of being seamlessly integrated with graphene in 
the fabrication of future nanoelectronics devices. One 
way to produce FGSs is by exfoliation from graphite ox- 
ide (GO), which can have different compositions with 
various oxidation levels depending on the synthesis pro- 
cesses and conditions. Currently, GO is of particular in- 
terest to scientists since chemical reduction of GO has 
been demonstrated as a promising solution-based route 
for mass production of graphene 

The electronic properties of the oxidized layers depend 
on the detailed chemical structure, which remains un- 
resolved for GO for more than a century with various 
structural models proposed in the literature 1^, 17. 18. 
[l9, 20]. Increased conductivity was observed during the 
reduction of an oxidized graphene sheet prepared from 
GO but the atomic origin of this behavior is still 
unknown. A recent high-resolution solid-state ^'^C-NMR 
measurement [2l| has confirmed the existence of C-OH 
(hydroxyl), C-O-C (epoxide), and sp^ C units on these 
layers. The data further indicate that a large fraction of 
sp^ C atoms are bonded directly to C atoms in the hy- 
droxyl and epoxide groups, and that a large fraction of C 
atoms in the hydroxyl and epoxide units are bonded to 
each other. In order to fully establish the atomic config- 
urations and related electronic properties in this impor- 
tant material, we report results from first-principles cal- 
culations that provide a clear picture of the energetically 
favorable building blocks and stable phases. Interest- 



ingly, various energy gap values are found for structures 
with different O concentrations, suggesting that the gap 
is highly tunable by varying the amount of hydroxyl and 
epoxide on the graphene sheet. 

Our theoretical study focuses on the following two key 
issues: (i) how these functional groups arrange them- 
selves on graphene, and (ii) how these arrangements 
affect the electronic properties of the graphene sheet. 
The calculations were carried out using the local-density 
approximation (LDA) within density-functional theory 
with a plane-wave basis set as implemented in the Vi- 
enna Ab-initio Simulation Package (VASP) Van- 
derbilt ultrasoft pseudopotentials |23i] are employed. All 
calculations are done with a plane-wave cutoff energy of 
500 eV. Results shown in Fig. 1 are obtained using a 
5x5 graphene unit cell with 5x5x1 fc-point sampling. 
The 5x5 unit cell contains 50 C atoms, providing a sep- 
aration of 12.3 A for the atomic combinations under in- 
vestigation. The vertical size of the supercell is 12 A, 
so the interaction between the layers is expected to be 
minimal. The size of the lateral unit cell is individually 
optimized for different coverages. The optimization of 
atomic positions proceeds until the change in energy is 
less than 1 x 10~^ eV per cell and the force on each atom 
is less than 0.02 eV/A. 

A single functional group of epoxide or hydroxyl on 
graphene can induce significant local distortion. With a 
new bond formed between C and O, the bonding char- 
acteristics of the connecting C atoms change from pla- 
nar sp^ to distorted sp'^ hybridization. The structure 
obtained in our calculation is in good agreement with 
that reported in previous theoretical studies [2J, |25[ . Of 
particular interest is the distribution of these functional 
groups on graphene. A recent atomic force microscope 
(AFM) measurement showed that the oxidized graphene 
sheets appear to have a thickness equal to integer mul- 
tiples of /i w 6.7 A [2^, indicating that epoxide and hy- 
droxyl are most likely to be present on both sides of the 
graphene sheet. Hence, all of our calculations have con- 
sidered possible two-sided configurations. 

To find out whether these functional groups prefer to 
aggregate or not, we have calculated the energy change 
associated with grouping these epoxide and hydroxyl 
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(b) 0+OH: E = -0.8 eV 
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(d) 0+20H: E = -2.6 eV 
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FIG. 1: (Color online) Atomic configurations and energies 
of various favorable combinations of epoxide and hydroxyl 
groups on the graphene surface. The energy shown is calcu- 
lated using a 5x5 unit cell and is with respect to the energy 
sum of well separated individual units on the surface. C and 
O atoms are represented by large spheres and H by small 
spheres. 



units in various ways, and the results calculated using a 
5x5 unit cell are shown in Fig. 1. The energies are refer- 
enced to the energy sum of separated individual units on 
the surface, and the most favorable, fully-relaxed atomic 
arrangements are shown in Fig. 1 for each combination. 
As one can see from the figure, the energies are lowered 
considerably when these epoxide and hydroxyl units are 
grouped together. Part of the reason is a cancelation 
of vertical structural distortion when these units can be 
located on both sides of the sheet. Significant energy 
gains are found if the OH units form 1,2-hydroxyl pairs 
(Fig. 1(c)) on opposite sides of the sheet; and the com- 
bination of 0-I-20H and 20+20H are also especially fa- 
vorable (Figs. 1(d) and (f)). In Figs. 1(d), (e), and (f), H 
points toward the neighboring O on the same side, yield- 
ing a configuration characterized by a hydrogen bond. 



The results in Fig. 1 conclude that these adsorbed units 
prefer to aggregate on the graphene surface. This is in full 
agreement with the experimental results inferred from 
the NMR data [H, [13, [2lj . 

We then explore ordered phases containing epoxide 
and hydroxyl groups incorporating the energetically fa- 
vorable building blocks in Fig. 1. Based on the energy re- 
sults, we only consider arrangements with the OH group 
appearing in 1,2-hydroxyl pairs. Each periodic phase can 
be specified by the relative amount of "free" sp^ C atoms 
(denoted by C* , corresponding to C atoms not bonded to 
O), epoxide (C2O), and the 1,2-hydroxyl pair [C2(OH)2 
with the connecting C atoms included]. The repre- 
sentative stoichiometry is C*_^_j^(C20)x[C2(OH)2]2^, or 
equivalently Ci+2:-(-j,02,(0H)2j,, with 0<a;<l, 0<y< 
1 , and < X + y < 1. The ordered phases we have calcu- 
lated are marked on a ternary diagram shown in Fig. 2(a), 
where the dashed lines indicate phases with the same ra- 
tio of epoxide versus the hydroxyl pair. For each phase 
we investigated, the structure was optimized by explor- 
ing various different local arrangements of the epoxide 
and hydroxyl groups. The lattice parameters and atomic 
coordinates are fully relaxed. The formation energy is 
defined in the usual way for a ternary system: 



AE[x,y] - S[Ci+,+yO,(OH)2,,] - (1 - 
^xE[C20] - yE[C2{OR)2] 



y)E[C*] 



(1) 



where E[Z] represents the energy of a periodic phase Z. 
These formation energies are shown in Fig. 2(b) for the 
"binary" phases and for phases along each dashed line in 
Fig. 2(a). 

For the fully oxidized phases in which all C atoms are 
bonded to O in either an epoxide or a 1,2-hydroxyl pair, 
low-energy "binary" phases with mixed epoxide and hy- 
droxyl compositions are found (top panel of Fig. 2(b)). 
These fully oxidized phases have a lattice expansion of 
the order of 2-3 %, and an LDA energy gap of 2.8 - 4.2 
eV. With a negative formation energy, these intermedi- 
ate phases are stable against separation into pure epoxide 
and pure hydroxyl phases. These stable phases include 
C602(OH)2, C40(OH)2, and C60(OH)4, with an epox- 
ide to hydroxyl ratio of 1:1, 1:2, and 1:4, respectively. 
As an example, the structure of C602(OH)2 is shown 
in Fig. 3(a), in which one can easily identify the follow- 
ing key features: the 1,2-hydroxyl pairs are connected to 
form a chain-like structure on both sides of the sheet in 
such a way that the interaction associated with hydro- 
gen bonds can be maximized; and O atoms are drawn 
to the remaining C atoms in the same hexagonal rings 
to form epoxides in the close proximity. The other two 
fully oxidized phases, C40(OH)2, and C60(OH)4 follow 
the same pattern for their atomic arrangements. These 
features turn out to be quite energetically favorable in 
constructing phases with intermediate compositions, as 
will be discussed below. The epoxide-only phase, C2O, 
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FIG. 2; (Color online) (a) Ternary diagram showing ordered phases on the graphene surface with different amounts of sp^ 
Carbon (C*), epoxide (C2O), and the 1,2-hydroxyl pair (C2(OE[)2). The phases investigated in this study are marked on the 
diagram, with dashed lines indicating those with the same relative amount of epoxide and hydroxyl pairs, (b) Formation energy 
as defined in Eq. (1) for different phases marked in the ternary diagram in (a), (c) Energy gap as a function of the overall 
oxygen-to-carbon ratio for different phases. 



is shown in Fig. 3(b), in which the O atoms follow the ar- 
rangements in Fig. 1(a) and stay in rows on opposite sides 
of the sheet [23| . The buckling is symmetrically compen- 
sated. For lower-concentration epoxide-only phases, the 
opening of the three-membered epoxide ring by breaking 
the C-C bond to release strain as proposed previously f24| 
turns out to be energetically favorable, together with the 
formation of O rows. 

Apart from the fully oxidized "binaries" shown in the 
top panel of Fig. 2(b), we are not able to find other or- 
dered phases with a negative formation energy. There- 
fore, the T=0 lowest-energy configuration of the oxidized 
graphene sheet is likely to be a combination of fully ox- 
idized regions and the clean graphene phase. This con- 
figuration has not been observed experimentally at fi- 
nite temperature. Possible reasons include the following: 
the entropy term could be important at finite tempera- 
ture; the oxidization process is a highly non-equilibrium 
one; and the covalent bonding in the epoxide and hy- 
droxyl units largely reduces their mobility on the surface. 
Therefore, domains of various intermediate phases may 
still be found in the sample under experimental condi- 
tions. The relative amount of epoxide and hydroxyl units 
on the graphene sheet depends on the sample preparation 
process and can vary over a wide range. Here we focus 
on phases along each dashed line in Fig. 2(a), in which 



the relative amount of epoxide and hydroxyl units on the 
surface is a constant. After exploring various atomic con- 
figurations, the lowest-energy periodic structure of these 
intermediate phases found in our calculation contains 
strips of epoxide and hydroxyl combinations with clean 
graphene ribbons in between. An example, the structure 
of C2402(OH)8, is shown in Fig. 3(c), which contains sep- 
arate regions of sp^ and sp^ carbon. The sp'^ strips con- 
sist of double hydroxyl chains and neighboring epoxides. 
The formation energies for the periodic phases with var- 
ious intermediate compositions constructed in this fash- 
ion are shown in Fig. 2(b). The strips interact weakly 
when separated by carbon ribbons, which explains why 
the formation energies of these phases fall on a straight 
line in Fig. 2(b) when the C* composition is larger than 
about 0.4. (Phases with different C* compositions have 
different strip separations.) With the formation energies 
falling on a straight line, many of these phases are ex- 
pected to coexist on the surface. It is energetically favor- 
able for these strips to coalesce, since the fully oxidized 
phase has a lower formation energy. However, this pro- 
cess may not be completed during the oxidation process. 

The current results indicate that the configurations 
investigated in previous density-functional calculations 
[2^ for oxidized graphene sheets containing both the 
epoxy and hydroxyl groups may not be energetically 




FIG. 3: (Color online) Atomic structures for selected phases, 
(a) Fully-oxidized phase of C602(OH)2. (b) Fully-oxidized 
epoxide-only phase C2O with oxygen rows on both sides of 
the plane, (c) C2402(OH)8 structure with hydroxyl-epoxide 
strips separated by sp^ carbon. C and O atoms are repre- 
sented by large spheres and H by small spheres and the dashed 
rectangles indicate the respective unit cells. The hydrogen 
bonds in the hydroxyl chains above the plane are indicated 
by dashed lines. 



favorable for a given composition. In our calculation, 
we find that the formation of hydroxyl chains result- 
ing from hydrogen-bond interaction between neighbor- 
ing 1,2-hydroxyl pairs greatly lowers the energy and 
that the epoxides are grouped next to these hydroxyl 
chains. The epoxide and hydroxyl units randomly de- 
posited on the surfaces are expected to arrange them- 
selves locally following these preferred patterns, giving 
rise to patches of sp^ carbon surrounded by fully oxidized 
epoxide-l- hydroxyl regions or vice versa. The possible ex- 
istence of pure graphene ribbons was proposed previously 
[2^ in order to explain the shift of the Raman peaks in 
GO and FGSs. The current study provides a strong sup- 
port for this picture based on extensive first-principles 
calculations. 

The finite region of sp^ carbon has a consequence in the 
electronic structure. Without knowing the exact atomic 
arrangements at various compositions, we use the results 
of the ordered structures investigated above to provide 
an estimate for this effect. The calculated LDA energy 
gaps associated with the ordered phases considered in 
Fig. 2(a) as a function of the overall 0/C ratio are shown 
in Fig. 2(c). The data-point symbols are the same as 



those in other parts of Fig. 2. An LDA energy gap of 
the size ranging from zero to 4 eV can be found for the 
range of O concentration we considered. The gap range 
includes both semiconducting and insulating phases, and 
the gap size is dictated by the width of the graphene 
ribbons in the ordered phases we studied. (It is well 
known that an LDA calculation underestimates the en- 
ergy gap, and the real gap value is expected to be larger.) 
A few vanishing band gap values are associated with arm- 
chair ribbons with 3n-|-2 rows of atoms (n is an integer) 
or zigzag ribbons with an even number of atomic rows. 
These may be considered as special cases. The overall 
band-gap results suggest a promising and practical way 
to tune the energy gap in FGSs by varying the degree of 
oxidization, which is feasible experimentally, and possi- 
bly the location of the oxidized regions. 

In summary, with first-principles calculations we have 
studied the structure and energetics of oxidation func- 
tional groups (epoxide and hydroxyl) on single-layer 
graphene, and the induced changes in the electronic prop- 
erties. We find that it is energetically favorable for the 
hydroxyl and epoxide groups to aggregate together and 
to form specific types of strips with sp^ carbon regions in 
between. An LDA band gap ranging from a few tenths of 
an eV to 4 eV can be obtained by changing the oxidation 
level and the location of the oxidized region, suggest- 
ing a great potential to tune the energy gap in graphene 
through controlled oxidation processes. 
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